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P-H ■ Abstract 

C ■ We discuss how an eigenvalue perturbation formula for transfer operators of dynamical sys- 

tems is related to exponential hitting time distributions and extreme value theory for processes 
generated by chaotic dynamical systems. We also list a number of piecewise expanding systems 
' to which this general theory applies and discuss the prospects to apply this theory to some classes 

^—■^ ' of piecewise hyperbolic systems. 



1 Introduction 



The first occurence of a rare event A e in a time discrete dynamical system T : M M can be 
described in terms of the first hitting time t £ (x) — inf{i ^ : T l x e A e } (x e M). In order to justify 
the adjective "rare" one should assume that n(A e ) is small, where fj, is a suitable reference probability 
measure on M. For many dynamical systems it is then known that the law of t £ under \x is nearly 
exponential with a parameter A e that is closely related to (J,(A e ), see [T8] for an excellent review. In 
this note we describe how asymptotic expressions (as e — > 0) for the exponent A £ of exponential hitting 
■ time statistics can be derived from spectral perturbation results for the Perron Frobenius operator of 

the dynamical system, provided the operator has a spectral gap on a suitable space of (generalized) 

o; 



functions. 



1.1 Perturbation theory for dynamical systems 

In 1999, C. Liverani and the author [IT] published a spectral perturbation theorem for linear operators 
1 that is particularly suited to deal with transfer operators of dynamical systems: it does not require 

the perturbation to be small in operator norm (which is rarely the case for transfer operators if the 
underlying dynamics arc perturbed). Of course there is a price to pay for this: the operators must 
satisfy a uniform inequality which is known as Lasota-Yorke inequality in dynamical systems or as 
Doeblin-Fortet inequality in the theory of Markov operators, and the conclusions are not as sharp as 
one might expect, namely the perturbed spectral quantities do not depend Lipschitz-continuously on 
the size of the perturbation. 

There are in fact examples showing that the latter drawback of this perturbation theorem is not 
just an artifact of the proof 8, Theorem 6.1 and Footnote 7]. Therefore, a perturbation result for 
the leading eigenvalue of transfer operators published by C. Liverani and the author [44] in 2009 
might seem contradictory at first sight: under rather general assumptions that are met in particular 
when the aforementioned spectral perturbation theorem applies, this result shows that the leading 
eigenvalue of a transfer operator is differcntiable as a function of the size of the perturbation at the 
unperturbed operator. This covers the case where the leading eigenvalue becomes smaller than 1 due 
to a small leak in the phase space of the system so that this eigenvalue determines the escape rate of 
the open system. 
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1.2 Escape rates in open dynamical systems 

The study of escape rates (also called decay rates) in open dynamical systems has different aspects. 
Starting with [50j [49] there are numerous publications highlighting the role these rates play for condi- 
tionally invariant measures, see [TSJ [33] for reviews and further perspectives. In [35] and [101 section 
9.6C] escape rates were characterized in terms of topological pressure, and [22] is a recent systematic 
account of this approach. 

Here, on the contrary, the focus is on precise approximations of the escape rate for small holes. 
This work was motivated by a result of Bunimovich and Yurchenko [TT1 Theorem 4.6.1] that I first 
learned about during a conference at the Erwin-Schrodinger Institute in 2008. It gave rise to the 
paper [44] . Since that paper was written, not only several other authors considered related problems 
on escape rates and intermittent phenomena [4J03I6, 22 24 s [28] [29] , but I also had the opportunity to 
learn about the close relations between escape rates, hitting and return time statistics and extremal 
indices for extreme value distribution^] [25[ \27\ 126] . In particular the extremal index turns out to 
be a quantity that appears also in the eigenvalue perturbation formula from |44] , and the sharp error 
term for exponential hitting times that was proved for mixing processes in [T] is also closely related 
to the same formula. 

1.3 Outline of this paper 

The goal of this note is to provide details for the interrelations described above in the particular case 
where the dynamics of open systems with small holes are described by a perturbed Perron Frobenius 
operator. In particular, 

— some details are given why the approximation formula for the leading eigenvalue is valid in the 
perturbative setting from [5T] (section^, 

— a slight extension of the arguments from [44] is given that provides estimates for the error term 
in the exponential approximation of the distribution of hitting times (section 13. 2p , 

— it is argued that exponential hitting time statistics (including a formula for the extremal index) 
follow from the formula for the leading eigenvalue (section [3.3[1 . 

— and a number of concrete settings are listed to which the perturbation formula applies (section 

2 A review of some perturbation results 

The paper [44] provides a first order approximation (as e — > 0) for the leading eigenvalues of the 
Perron- Frobenius operator of an open dynamical system with a hole of "size" e. The main result is 
formulated as a perturbation theorem for the leading eigenvalue of a family (P e ) eg e of linear operators 
on a Banach space (V, ||.||), where E is a set of parameters (equipped with some topology) containing a 
parameter called such that Pq describes the dynamics without a hole. This theorem applies whenever 
the following assumptions ^ - ([6]) can be verified: 

There are A 6 € (D, <p e 6 V, a bounded linear functional v £ : V — > <D, and linear operators 
Q £ : V — > V such that 

A~ P e = tpe ® v e + Qe (where we assume A = 1), (1) 
P £ (cp £ ) = \ £ <p £ , u £ P £ = X £ u £ , Q £ (cp £ ) = 0, v £ Q £ = 0, (2) 

OO 

V sup || Q^ 1 1| =:Ci <co, (3) 



1 I want to mention in particular the talk by J. Freitas "Laws of rare events for chaotic dynamical systems" and 
stimulating conversations with L. Bunimovich and M. Demers at the conference "Large deviations in dynamical systems" 
at the CIRM in June 2011. 
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Observe that assumptions ([TJ and ([2]) imply v £ {tp £ ) = 1 f° r a ll £@ The fourth assumption relates the 
"size" of tp £ to that of ip?,: 

3C*2 > OVe G E : vo((p £ ) — 1 and \\ip s \\ < C 2 < oo. (4) 

Now let 

A £ :=^ (P -P e )(^o), 

where fo(Po — P e ) denotes the composition vq o (P — P) : V — >• (D. Then the last two assumptions 
are 

]im||i*(i%-P e )||=0, (5) 

\\MPo-Pe)\\ •||(Po--Pe)^o|] < const -|A £ | (6) 
Under these assumptions the main result of |44j is the formula 

1 - A £ = A £ 0(1 + o(l)) in the limit e -» (7) 

where 6> is a constant that takes care of short time correlations, see section [3~T1 for more details. 

2.1 Verifying assumptions ([I]) — ([4]) 

One way to verify assumptions (UJ - (J3J that was already taken in [33] and [23] > is to use the spectral 
perturbation theorem from (31]. It holds (not only) under the following assumptions: 

i. There are constants a £ (0, 1), D > and a second norm \.\ w ^ |.| on V (it is enough that this 
is a seminorm) such that: 

The residual spectrum of each P £ is contained in {z £ (D : \z\ ^ a}. (8) 

Ve G £ V/ G V Vn G N : |P £ n /|» < ^ \f\ w (9) 
Ve 6 P V/ G V Vn G N : ||P e "/|| < Pa n ||/|| + D|/|„ (10) 

ii. There is a monotone upper semi-continuous function ir : E — > [0, oo) such that 

lim tt £ = and V/eVVee£: |P £ / - P /U < 7r e ||/||. (11) 

An equivalent way to rewrite the last condition is: 

Ve£E: |||P £ -Po|H7r £ 
where |||P||| := sup{\Rf\ w : f £ V, ||/|| 1} for any linear operator R : V -> V. 
The main result of [31] is the following theorem: 

Theorem 1 Suppose that (P £ ) £< ze is a family of linear operators on V satisfying assumptions |3|)- 
111)) . Fix 5 > and r £ (a, 1), and let p := log ^/ log — G (0, 1). Then there are constants a = a(r) > 
0, b = b(5, r) > 0, c = c(5, r) > and d — d(5, r) > and a neighborhood Es, r of in E such that for 
e £ Es^r and z £ (D with \z\ > r and dist(z, cr(P )) > S, 

\\(z - P £ )-7ll < a||/|| + b\f\ w for all f G V (12) 

and 

III (z P.)" 1 -(z- Pa)" 1 !!! < n§ ■ (c\\(z - Pom + d\\(z Po)- 1 !! 2 ) . (13) 

2 Assumption J3} is clearly equivalent to the seemingly stronger assumption that sup £gB decays exponentially 

as n — > oo. The form of assumption 10 chosen here underlines however, that formula {TJ can be proved as well under 
slightly different assumptions that allow in particular to use a weaker operator norm in (0 which is not submultiplicative, 
see gU Remark 2.2], 
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Remark 1 In [41] it is assumed that E = [0, oo), but as the parameter e enters the proofs only via 
the numbers 7r e , the same proof applies to a more general set E of parameters. Similarly, in |44j it 
is assumed that E is a closed subset of It, but again the parameter e enters the estimates only via 
the derived quantities ry £ and A e defined in section [3.11 below, so that also this result is valid for more 
general sets of parameters. 

Remark 2 If assumptions ([9]) and (TTO]) are met and if the closed unit ball of (V, ||.||) is \.\ w -compact, 
then <r(P e ) \ {z G (D : \z\ ^ a} consists only of isolated eigenvalues of finite multiplicity so that 
assumption ([8]) is automatically satisfied. This is guaranteed by a strengthened version of the Ionescu- 
Tulcea/Marincscu theorem (see e.g. [55]). 

There are numerous situations where this remark applies: 

— V is a space of Holder continuous functions on a compact metric space and \ .\ w is the supremum 
norm. If the underlying space is not compact, weighted norms may be used; reference |33j 
provides many examples. In any case this setting is adapted to the study of Markov operators 
acting on observables and to Ruelle transfer operators, see Ruelle's monograph [51, Proposition 
5.24] and [20] for a recent contribution. As this setting is less useful for studying open systems, 
no further details are provided here. 

— V is the space of functions of bounded variation on a regular domain of finite volume in It d and 
\.\ w is the LLeb" semmorm - This setting is used for transfer operators of piecewise expanding 
maps, starting with [37l [38l [36l [52] and being continued in numerous papers, see [7| for an 
overview until 2000, [21] for references until 2008 and [54] for a very recent contribution to 
this direction. In practically all variants of this setting the passage to an open system can be 
described as a small perturbation of the Perron-Frobenius operator in the sense of the |||.|||-norm. 

— More generally, V is a space of distributions on a finite-dimensional smooth manifold, equipped 
with a suitable norm, and \.\ w is another weaker seminorm on V. Starting with JT5j, such 
spaces have been applied so far to various hyperbolic or piecewise hyperbolic systems, see e.g. 

[321 queue ma. 

— V is a space of distributions on [0, l] z which contains in particular the space of signed measures 
with absolutely continuous marginals of (uniformly) bounded variation (see [121 Ell S3] for 
details). This setting was used to study weakly coupled lattices of piecewise expanding interval 
maps. Inequality ([TO")) was stated explicitly in Qj2j Lemma 3.4]. 

Remark 3 The conclusions of Theorem Q] remain valid, if assumption ([8]) is replaced by the following: 
If S > and r — a > arc sufficiently small, then the set 

{z G (D : \z\ > r, dist(z, a(Po) > 5} is connected., 

see [41] Remark 6] . 

Now we specialize to the case where the following operator mixing condition holds: 

Po = Vo ® Vq + Qo as m © where Q has spectral radius strictly smaller than 1. (14) 

Observe that this implies the ([8]) for some a G (0, 1). Therefore, assumptions ([9]) - ([PI]) and flH} imply 
the validity of ([TJ) - J3]). This follows easily from Theorem [l] see [41] Corollaries 1 and 2]. It remains 
to show that also ([!]) follows: Observe first that by [2U Corollary 1], 

\(p £ - tp ■ v (tp e )\ w = \ip e ® v e (ip e ) - ip ® v (tp e )\ w Knr§ ■ \\<p e \\ < K 1 K 2 ttP ■ \cp s \ w 

for suitable K\ , K 2 > and p from Theorem [TJ Hence 

\tp e \ w ■ (1 - K 1 K 2 ir p e ) ^ U • ko(^ e )| 



4 



so that \ip e \ w < 2 • l^oliu • l^oCVe)! when e G E is sufficiently close to 0. Invoking [4TJ Corollary 1] once 
more, this yields 

||v?e|| ^ 2i^2 • • l^oCVe)! for £ £ E sufficiently close to 0. 

Therefore one can normalize tp e (and in consequence v £ as well) such that Vo(<p s ) — 1, and (j4|) holds 
for these normalized ip £ with a constant C2 = 2K2\(po\ w - 

2.2 Verifying assumptions ([5]) and (]6]) 

In the remainder of this paper we restrict attention to the following dynamical setting: 

REPFO (Rare event Perron- Frobenius operators) V is a space of functions or distributions on 
a metric space M that contains the constant functions, and P$ is the Perron-Frobenius operator 
of a dynamical system T : M — > M w.r.t. the reference measure vq . It is mixing in the sense of 
(fT4|) . The family (A e : e G E) is a family of subsets of M which are sufficiently regular such that 



i. the operators P e , P E f — Po(l-M\A e f)i satisfy assumptions ((8} - (fTTj) and (fl4"|) . 
ii- lAe-f€V for all / G V, and 

hi. |z/ (l^ e /)| • ||U e <A)|| < const ||/|| • |i/ (l^Vo)| for all / G V. 

Remark 4 If V is a space of functions such that sup |/| ^ ||/|| for all / G V, if ^0 is a measure on 
M, and if inf^gy a e (fo(x) > 0, then REPFO. hi is satisfied. 

Remark 5 As vqPq = fo, condition REPFO. hi implies (in fact, is nearly equivalent to) 

\\MPo - Pe)\\ ■ \\(Po - Pe)(^o)\\ < \MP0~Pe)Ml 

Let (Uo be the probability measure on M defined by /j,q(A) = J A ipodi'Q. (It is the only T- invariant 
probability measure absolutely continuous w.r.t. Uq.) Then 

HoiA) = MPo ~ Pe)(<P0) < COnst \\M P ~ P e)l 

Proposition 1 In a REPFO setting, assumptions (QP~([IP are satisfied. 

Proof In the preceding section it was argued that assumptions (JTJ) — d3J* are satisfied. Next observe 
that for f £ V with ||/|| < 1, 

\\fo\\ ■ MW)| = IK^o ®^o)(UJ)|| 

= lim ||P "((P -P E )(/))|| by© 

< D 3 \(P - P e )(f)\ w by (HUD 

< TTeDsll/H <7r e U 3 . by CO 

Therefore, 

WMPo - Pe)\\ = Sup MW)| *S *eT^ 

ll/IKi II || 

and 

WMPo - p e)\\ • \\(Po - PeWoW < const sup \p (1 A J)\ • ||U^o|| by (USD 

ll/IKi 

s$ const |i/o(1a s <Po)| by hi. 

= \v {P - P s )(tp )\. 

This proves © and ©. □ 
In section [4l we discuss examples of REPFO settings. 
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3 Rare events, exponential hitting time distributions, and ex- 
treme value statistics 



Recall that we restrict to the REPFO setting in the remainder of this paper. 

Our approach to estimate first hitting time distributions of rare events proceeds along the following 
well known lines: denote the first hitting time to A £ by r e , i.e. let T e (x) = inf{i ^ : T l x £ A e } for 
x £ M. We are interested in estimating vq{t £ ^ n} and fio{ T e 3 s n }- 

If tp € V is such that ipi/Q is a probability measure on M, then 



« n— 1 
Jm \i=0 



<y3 dj/ 



= A; 1 / tp e ® i/ B (y>) di/„ + A™ / Q^d^b 
Jm Jm 

= \y e (<p) + o(\«\\Q n j-y\\) (is) 

Here the constant implied by ,,0" is uniform in n and e. Applied with ip = 1 or ip = </?o one obtains 
expressions for fo{T e n} and ^o{ r e 3 s n }j respectively. To make use of this identity, we need some 
control of X e as a function of e. 

3.1 The perturbation formula for the leading eigenvalue 

The main result of 44j asserts that A e = exp (— fio(A e )6(l + o(l))) with a constant £ (0, 1] deter- 
mined by the dynamics and a o(l)-term for s — >■ 0. 

For the further discussion we recall a more precise intermediate identity from [44] . Let 

% := |MP -P e )||= sup y fdvo-.feV, ll/H 

A £ := vo(Po - P e )(<Po) = / <^o = Vo(A e ), 



£sup||Q*|| = 0((l-7) JV ), (16) 



where 7 > is some lower bound for the spectral gaps of the P e . (It is strictly positive because of 
Theorem[TJ see also [3U Corollary 2].) Then we have, for each e £ E with A £ ^ and for each N £ N, 

= (^, E + - 7 )")) (1 + 0(N % )) (17) 

with constants independent of N and e [321 eq. (6.2)jf|, where 

»*.-!- -w. ^.- ^w-^ww-ftxw)) . (18) 

fe=0 e 



3 If A e = 0, the numbers q k e are not well defined. Therefore the argument, given in 1441 after eq.(6.2)], that 
A e = Aq in this case, is wrong. The claim is nevertheless true, because assumption ifjJJ implies that vq(Po — P s ) = 
or [Pa - P e ){<p ) = 0. In both cases, (A /A £ )™ = (A / 'X e ) n M<Po) = u (X7 1 Pe) n (fo) = M<Pe)M<Po) + o(l) + when 
n — > 00 in view of assumptions {1} - (3), so that A e = Aq in view of Jit . 
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For REPFO-type P e we are looking at here this simplifies, in case fio(A e ) > 0, to 

i f , XJ no(A e nT- 1 A c £ nT- 2 A c £ n...nT- k A c e nT- ( - k+1 '>A e ) , . 

Observe that J2T=o 9 fe .£ = 1 ^ or cacn f G E by Kac's recurrence theorem. On the other hand, 
assuming that 

oo 

<7fc := lim <j% £ exists for all k and denoting 6* := 1 — \ q/., (20) 

fc=0 

this latter sum can attain any value between and 1. In any case, if these limits exist, and if wc 
denote 9 N = 1 - J2k=o 9k, we g et 

As this is true for every TV, we conclude, as in [44], that 

1 - A °° 

SAUSS)-^ 1 -^ (21) 

i.e. 

A £ = 1 - 9vo(A e ) + o(/i (A s )) = exp(-^ (A £ ) + o(fx (A e ))). (22) 

In the next subsections we will draw some finer conclusions on the (non)-hitting probabilities vq{t £ 
n} and /xo{r e ^ n} working directly with equation (|I7|) . 

Remark 6 For the map Tx = 2xmodl and vq the Lebesgue measure on [0,1], formula (|2"Tj) was 
derived by combinatorial methods in [17l Theorem 4.6.1]. In the context of Gibbs measures with a 
Holder continuous potential this formula was derived also in [24], using spectral theoretic methods 
and ideas on exponential first return times from [34| . 

Remark 7 If the map T is continuous, if the holes A e shrink to a single point x G M and if fj,o(A e ) > 
for all e 7^ 0, then 

i. qk = lim e _>o ^^^"o ) — if x is a periodic point of minimal period k + 1, 

ii. qk = otherwise. 

Therefore, 9 < 1 is possible only if x is a periodic point. If x is a hyperbolic periodic point, then 
9 < 1, indeed. 

Remark 8 Let T : [0, 1] — > [0, 1] be a piecewise monotone interval map of Rychlik type (which always 
fits the REPFO setting, see section l4~Tj) . Denote by i^o the Lebesgue measure on [0,1]. Assume 
furthermore that V is a finite subset of (0,1), that A £ is the e-neighbourhood of V, and that the 
invariant density ipa is continuous at each x S V. 

Denote by II the set of all pairs (x, k) e V x N such that T k+1 x <E and T j x g" V for j = 1, . . . , k. 
It is easy to prove that 

' = 1 " ( ,5nl^)T' 
It follows that < 1 if and only if the set II is nonempty. See also the related examples in 03] . 

Remark 9 If the holes shrink to a lower dimensional submanifold which is invariant under the dy- 
namics, then typically qo > (and hence 9 < 1), see [U] section 4] for a particular example. 
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3.2 A sharp error term for exponential hitting time distributions 



Abadi [l studied the convergence of hitting time distributions for hitting cylinder sets in a shift space 
over a finite alphabet A under 0-mixing stationary laws. He proved: Let /iq be a 0-mixing stationary 
measure on A N with J2 n < M n ) < 00 • Then there are constants C > and 0<Si<l<S2<oo 
such that for all LeN, cylinder sets A = [o,q . . . 0£_i] and t ^ 1, there exists £a G [Si, S2] with the 
following property: for the first hitting time ta of A and all t > 0, 



Mo i t a > 



C5(A)(tVl)e- 



(23) 



where 5(A) = mm L ^ N ^ 1/llo{A) (Nfi (A) + (f>(N)). 

An analogous result follows in the dynamical REPFO setting. One should observe, however, that 
Abadi assumes only a summable mixing rate whereas the REPFO setting always implies exponential 
decay of correlations for "regular" observables (i.e. for observables in V). 

Proposition 2 Suppose that in a REPFO setting fio(A £ ) > for all e £ E and that the number 8 
defined in \21)) is strictly positive. Then there is a constant C > such that for all e G E sufficiently 
close to there exists £ £ > with the following property: for the first hitting time t £ of A e and all 
t > 0: 

t 



Mo S T e ^ 



s£ C6 E (t V l)e" 



where S £ — inhii^N (Nr] £ + ^jv) 



&Mo(^s). 

0(r\ £ logrfe). The numbers £ e satisfy lim e _ ! .o £ e 



(24) 



Proof Given e > 0, fix N 6 N such that iV?7 e + O(kn) becomes minimal, and denote this quantity 
by S e . Here 0(n N ) = 0((1 - j) N ) is the term from equation (02). Then TV = 0(logr; £ ) and 

5 e = 0(lh logjfc). Let £ E = 6>jV,e + Then 



A e = 1 - Mo(A e )£e (1 + 0{Nn £ )) = cxp[-fi a (A £ % (1 + 0(5 £ ))} 

as Mo(^e) < |[ VO [| • Ve- 

On the other hand, X e = 1 - /x o OA e )(0 + o(l)) when e by ([2l]). Therefore 



(25) 



(e-e e )(i + o(« e )) = eo(« e ) + o(i) 







so that lim e ^o £e = > 0. In particular £ £ > if e is sufficiently close to 0. 

Now let n = \ t/(^ £ ^(A £ ))\ and assume for the moment that n > 1 losv 
and (fl"6"|) . there is some q S [0, 1] such that 



£eMo(^e 



= A^K(^o) + 0(\\Q?\\)} from(H3 
= exp[-t(l + 0(4))] • (1 + Ofafe)). 



iog(i- 7 ) i- Then ' in view of 0B 

(26) 



(Here we also used v £ {fo) = 1 + 0(r] s ), see [44j Lemma 6.1].) Therefore, 

' < C(tVl)e^ E 



(27) 



(For t S e ^ c, where c > a suitable small constant, one can estimate by C e t S £ , in the other case 
by c- x Cte-H £ ) 



For n ^ 



log Si 
log(l-7) 



we have the following trivial estimates: 



1 > Mo {r £ > - 1 — X > MoK > n + 1} > 1 - Y)/io(T-%) = 1 - (n+ l)Mo(A e ) ^ 1 - 0(4) 



and 



1 > e~* > 1 - * > 1 - (n + l%(Mo(Ae) > 1 - 0(<J e ). 



Hence the l.h.s. in (f2"7) is bounded by 0(<5 e ), and as e 4 ^ 1 — 0(<5 e ), this proves (|2"T1) also for 



log (1-7) I 



□ 
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Remark 10 The problem of convergence of the distributions of rescaled hitting times to an exponen- 
tial without an error term as in (|24[) was addressed in many papers and under various assumptions, 
mostly for holes A e that are not necessarily cylinder sets. Galves and Schmitt [30] started these 
investigations under an exponential ^-mixing assumption. Hirata obtained a slightly weaker form of 
such results for Gibbs states on subshifts of finite type and for Axiom A diffeomorphisms [53] and for 
system with a finite generating tp- mixing partition |35j as an important step in proving Poisson laws 
for these systems. For a-mixing prcesses over a finite or countable alphabet, Abadi and Saussol [3] 
proved convergence of the distributions of the rescaled hitting times to an exponential without any 
further assumption on the mixing rate. Abadi [2] discussed in some detail the nature of the scaling 
factors 

3.3 Extreme value statistics and the extremal index 

In a series of papers, Freitas, Freitas and Todd ES1 ESI EI] studied extreme value statistics for obser- 
vations generated by chaotic dynamical systems and related them to first hitting time statistics. The 
starting point of their approach are mixing conditions in the tradition of Leadbetter [351 ES] that are 
adapted to the problem of studying extremal events. They show that these conditions are in particular 
satisfied if the underlying system is a (topologically mixing) one-dimensional map of Rychlik type |52j 
and vq and [1q are the associated conformal resp. invariant measure. Here we show how some of their 
results can be derived (and slightly extended) from our formula (|2ip . 

We begin recalling the approach to extreme value laws taken in j25j HSJ [27] . Suppose that X : 
M — > R is a continuous function. Let z max := sup{z £l:!/o(I)z)>0} and assume that for some 
za < W all sets V z :— {X > z} for z € [zq, z max ] satisfy the assumptions on the sets A £ from above. 
Suppose for simplicity that fo(V ZjxlM ) = 0. Fix t > and determine real numbers z n G [zo,z max ] 
such that lim^oo n/j, {X > z n } = t. Consider the process X n = X o T n on the probability space 
(M, B, v ). Then max(X (u;), . . . , A„„i(w)) z n if and only if T 1 uj e M \ V Zn for i = 0, . . . , n - 1. 
Denote e = e n := z max — z n . Then equations (| 15[) and (|22p imply the following extreme value law 

z/ {max(Xo,...,X n _i) < z n } = v q {lu e M : T 1 lj <£ V Zn for i = 0, . . . , n - 1} 

= Mn^n} (28) 
= exp(-i0(l +o(l))) 

in the limit n — > oo, where also the fact that lim £ ^o t/ e(l) = ^o(l) = 1 (see [41]) is used. In this 
context the constant 8 defined in (|20|1 is called the extremal index. It follows from Remark [7] that in 
case T4 max = {x} is a singleton, 9 < 1 is possible only if x is a periodic point. In view of the results 
collected in section 0] this answers a question by Freitas, Freitas and Todd from the introduction of 

ESJ. 

4 A review of open systems to which the eigenvalue pertur- 
bation formula applies 

In this section we discuss some classes of open systems of type REPFO. The results of section [3J apply 
to these systems. 

4.1 Piecewise monotone interval maps of Rychlik type 

Rychlik |52j provided a very elegant and useful scheme to treat piecewise monotone one-dimensional 
maps T : [0, 1] — > [0, 1]. The ingredients of his approach are a non-atomic probability measure vq on 
[0, 1] with full support and a function g : [0, 1] — > [0, +oo) of bounded variation. It is assumed that 

a) there is an open set U C [0, 1] with vq{U) = 1 such that for each connected component J of U, the 
restriction Tjj is a homeomorphism from J onto TJ (denote the collection of open components of 
f/by/3), 

b) g\ s = where S := [0, 1] \ U, 



c) ■& := lim, woc llfifnllco" < 1 where g n = g ■ g o T ■ . . . ■ g o T n l , and 

d) V(j[Pf) = vo(f) for each bounded measurable / : [0, 1] — > 1R where 

Pf(x)= J2 a(y)f(y)- 

yeT- x x 

Denote by BV the space of functions of bounded variation on [0, 1] and define the seminorm \f\ w = 
J \f\dv and the norm ||/|| = Var(/) + \f\ w on BV. Then assumption d) implies (j9]) for e = 0, and 
Proposition 1 of [55] shows that 

Va G (ti, 1)3C > 0V/ £ BVVn e¥S : Var(P n (/l B )) ^ On 11 Var(/) + C f \f\dvo (29) 

Be/3™ J 

where /3 n = /3 V T^fi V ... V T^" -1 )/?. This is (10]) for e = 0. Finally, © for e = follows from [5j] 
Theorem 1]. 

Let z £ (0, 1) and consider intervals A e = (z — e,z + e) C (0, 1). If we replace S by S U {z}, all 
of the above assumptions a) - d) remain valid. So w.l.o.g. z £ S. The rare event PF operators P £ 
satisfy 

\PefU = f \f\dva < f \f\dvo = |/U , (30) 

• / n;ir-*([o,i}\^) ■/ 

and, observing that £> n Dfeo r _fc ([0, 1} \ A e ) is a subinterval (possibly empty) of B for each B £ /?" , 
also 

Var(F;7) < J! Var(F £ "(/l B )) £ Var(P"(/l B )) < Ck" Var(/) + C\f\ w . (31) 

Be/3" 

The uniform estimates ([30]) and (|3T]) show that © - (| 10[) are satisfied for all e > 0. We turn to (fTT]) : 

|P E / - Pf\w < / |/|^0 < ^o(^)||/||oo < ^0(A £ )||/||. 

So r e = ^o(^-e) \ as e — ► 0, and this is (ITT]) . Thus assumption REPFO.i is verified. 

Assume now that T is mixing in the sense of (|14p with unique invariant density ipQ. As tpo is of 
bounded variation, the following extra assumption is not very restrictive: 

Positivity assumption: S := lim e _>.o max {inf (z- e ,z) Po, ' m ^(z,z+e) Po} > 0. 

Indeed, if T has only finitely many monotone branches, Theorem 8.2.3 in |15| shows that this assump- 
tion holds if z belongs to the support of the invariant measure fo vq. 
Under the positivity assumption, 

||WII= / \f\duo + Var(l A ,/)<|/| w +2|[/[[ 00 + Var(/)<3||/|| 

J A e 

and 

MW)| ■ WU.foW < MAe)\\f\\oo • 3 ll^oll < 4<5 1 ||^o II • vo(1a.<Po) • 11/11 
for all suficiently small e > 0. These two estimates yield REPFO.ii and iii. 

Remark 11 The class of Rychlik maps contains the classes of maps considered in [JJ and [IS]. Freitas 
et al. [25] proved the validity of an extreme value law for Rychlik maps along the classical 
probabilistic lines initiated by Leadbetter (45j [46] . Exponential statistics for Rychlik maps around 
fio-a,.e. point z were already derived in [16] . 

See also Remark |5J for an explicit formula for 9. 
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4.2 Piecewise expanding maps in higher dimensions 

For piecewise expanding maps in higher dimensions, the classical BV-space is the first candidate for 
a suitable function space V (3TJ [37l [54] . For the present purpose, however, the quasi Holder spaces of 
Blank jI3J and Saussol [S3] are more suited, because, in contrast to classical BV-spaces, these spaces 
are algebras and their norms control the L°°-norm [53j Proposition 3.4]. 

Without going into details, we consider a piecewise C 1+a map on a compact set SI C WL N for 
which expansion wins over the effects of discontinuities in the sense of ;53, condition (PE5)]. The 
corresponding PF operator P satisfies (JHJ - (fTO]) for e — in view of [53j Lemma 4.1 and Theorem 
5.1]. Fix a point z in the interior of the domain Ui of one of the branches of T, and denote by A e the 
e-neighborhood of z (w.r.t. some £ p -norm in IR ). Then the corresponding rare event PF operators 
P e satisfy the same estimates that P satisfies if the domain Ui is replaced by U \ A e (for sufficiently 
small e). This modification of U increases the constant G in (PE5) of |53jj by at most const 
and it follows from the explicit expressions for the constants in (|10j) that are given in the proof of 
Lemma 4.1 in [53], that ([HJ - (fTOj) are indeed valid for all sufficiently small e > 0. Assumption (fTTj) is 
satisfied because ||/||oo ^ const ||/|| where ||.|| is Saussol's quasi Holder norm. This is REPFO.i. As 
\\l Ae \\ is of the order of e^" 1 , REPFO.ii follows from ||U S /|| < const ||1 A J • ||/||, see [5j Proposition 
3.4]. Under the additional positivity assumption S := lim^o infoi e ipo > 0, condition REPFO.iii 
follows immediately from ||/||oo ^ const ||/||. 

In higher dimensional systems, holes that do not shrink to a point may occur quite naturally, for 
example when A £ is defined as the event that a state vector has two nearly (i.e. up to e) identical 
coordinates. A simple example of this type, namely a direct product of a piecewise expanding interval 
map T with itself, is discussed in |19j . where exponential hitting time statistics are proved (as a 
corollary to a Poissonian limit law) and where a (complicated) formula for the exponent of the limit 
distribution is derived also in the case where only some power T m of the map is uniformly expanding. 
In |44[ section 4] this formula is rederived from spectral perturbation theory when T itself is expanding. 
Indeed, for m > 1, too, this formula can be derived from (|2ip by a straightforward, though lengthy 
calculation. In the same reference also a modification of the product system by some weak coupling 
is discussed. 

4.3 Gibbs measures on subshifts of finite type 

Ferguson and Pollicott [24] study escape rates for Gibbs measures with Holder continuous potentials 
on subshifts of finite type. They obtain formula (|21|) for suitable holes (not just cylinder sets) in 
this context and use it to derive from that the corresponding formula for the escape from conformal 
repellors relative to equilibrium states, and also a related formula for the Hausdorff dimension of the 
set of points on the repellor that never enter the hole. For deriving formula (|2Tj) they produce a 
function space and two norms fitting the REPFO setting. 

4.4 Topological entropy of perturbations of shifts of finite type 

If a subshift of finite type over a finite alphabet is perturbed by deleting one or several blocks and if 
these blocks have small mass under the measure (1q of maximal entropy of the unperturbed subshift, 
the topological entropy of the systems drops by a small amount. As the topological entropy is the 
logarithm of the leading eigenvalue of the corresponding transfer operator with constant weight 1, 
the amount by which the entropy drops down is controlled by formula (|21[) when it is applied to the 
operators normalized such that the unperturbed operator has leading eigenvalue 1. 

Lind [47| Theorem 3] provides upper and lower bounds on the entropy difference when a single 
(long) block is discarded. These estimates can be made more precise using formula (|2"TT) . Details of 
this in two slightly different settings can be found in [44] section 5.2] and in [24] Corollary 5.4]. 

4 The constant there is denoted G(e,£o); note that our e (size of the hole) and the e in G(e,£q) are unrelated. 
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4.5 Further candidates 



The investigation of statistical properties of dynamical systems via the spectral properties of their 
associated Perron Frobenius operators is currently a very active area of research. Some of the available 
results apply to piecewise smooth systems and there is good hope that these systems are of REPFO 
type. Here we discuss briefly the situation for piecewise hyperbolic maps and for coupled map lattices 
of piecewise expanding ID maps. 

Piecewise hyperbolic maps: Demers and Liverani [21] studied Perron-Frobenius operators for a 
broad class of piecewise hyperbolic maps on two-dimensional compact Ricmannian manifolds. They 
provide a normed space (£>, ||.||) of distributions and prove in [211 Lemmas 3.5 and 6.5, Proposition 
6.6] that assumptions © - (jTTJ) hold for many reasonably regular holes A £ . In particular this is true 
for disks of radius e when s is sufficiently small (so that the curvature of dA E is large, see [2TJ Remark 
2.16]). Hence REPFO. i is satisfied. A first step towards a proof of REPFO. hi is the observation 
that |fo(lA £ /)| ^ £l+a ||/l|s ^ ll/H for all / € V according to the definition of the stable norm 
||.|| s and the strong norm ||.|| in (21[ eq. (2.3)and (2.5)]. Therefore, it mainly remains to show that 
e 1+a \\l As ip \\^v (l Ae ip ). 

Baladi and Gouezel [HI [ID] studied the Perron Frobenius operator of piecewise hyperbolic systems 
in any dimension on certain classes of Triebel spaces which they chose such that multiplication by a 
regular set, e.g. a convex set, acts as a bounded linear operator. This indicates that also their systems 
might fit into the REPFO setting. 

Coupled map lattice of piecewise expanding interval maps: Coupled map lattices (over 
Wi d ) of piecewise expanding interval maps were studied in [421 143] using a spectral theoretic approach. 
Here V is a space of distributions on [0, l] z that contains in particular all probability measures whose 
finite dimensional marginals have densities of bounded variation with uniform bounds on the variation 
of all these marginal densities in all coordinate directions. While [42] treats only coupling mechanisms 
which are uniformly weak on the whole phase space, the short note [43] deals with cases where strong 
(and discontinuous) coupling effects are allowed on small parts of the phase space. The estimates used 
in this context should be helpful also in attempting to prove REPFO. i-iii. 
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